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i^BSTRACT 


The analysis of composite plate is an increasing 
interest in aerospace structxires. Various methods avail- 
able in the literature have limitations interms of compu- 
tational efficiency* It is considered as an important 
criteria in optimal design* Keeping in view* a study has 
been made to apply successive iteration method in the an? 
sis of laminated plates. 

By Levy’s solution procedure* the governing part^ 
differential equations of the plate are converted to ord: 
differential equations* The iteration technique is then 
applied to eigen value plate problems# Both couploi and 
decoupled antisymmetric laminates have been considered i 
the present study. 

In decoupled laminates* biaxial buckling problem 
are investigated in cases with x edges simply supported : 
and y edges in simply supported* simply gupported-cl amp^ 
clamped* cl araped-f ree# simply supported-free and free 
boui^ary conditions, free vibration problems are also 
investigated in cases with x edges simply supported andi 
ybedges simply supported* clamped* and free boundary 
conditions# Conputations were carried out on Glass-epoj 
and carbon-epoxy homogeneous laminates for different aaj 



ratios, fibre orientations and biaskling load ratios. The 
eigen values# buckling load or natural frequency obtained 
using this technique agrees closely with the earlier 
methods* The technique consumes less# computational time 
and requires less efforts* The technique is independent < 
initial starting ftmctions* Optimal design of thin larain^ 
can be made readily# applying this technique with decoupl 
condition* 

la coupled laminates uniaxial buckling problem is 
investigated in case with x edges simply supported and y < 
clamped condition* Various computational model schemes w 
three and two dependent coupling variables failed to give 
correct result# eventhough convergence in eigen value was 
achieved* This may be due to sensitiveness of the coeffii 
of the derivatives appearing in the governing equations o: 
laminate* 
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GHAPTER-I 


lOTRQDUCTIQN AHD LITERATURE SURVEY 

1,1 

Fibre reinforced composite plates have been widely 
tised in the primary structures of aircrafts and spacecrafts 
because of their high specific strength and stiffnesses. 

The anisotropic properties can be tailored to the structiirat 
requirements by a judicious selection of lamination stacking 
sequence# fibre orientations# number of plies aid thickness 
of individual lamina. This flexibility in design leads to 
overall weight redijction. The other advantages such as 
durability# ease of repair# high energy absorption and seam- 
less construction offer an attractive alternative to metallic 
aircraft structural components. 

Composites exhibit shear extension and bending- 
stretching coupling# unlike isotropic materials. Couplings 
are undesir^le as it causes warping of laminates and 
increases ■ the complexity in analysis. Their presence reduces 
strength and stiffnesses of laminates considerably. It is 
possible to eliminate the effect of coupling by suitably 
manipudating thicknesses# fibre orientations and stacking 
^sequence. For optimal composite plate design decoupled • 
laminates may have to be used. An efficient technique is : 
desir^le for optimal design of composite plates# which tend 



2 


to decouple at the optimum point ^ A io tailed literature 

survey has bean made in the ted'iniques available for the 
analysis of composite plates. 

1.2 LITERATURE SURVEY 

T;'i.mcs.hdako and Krieger LlU studied the buckling of 
isotropic homogeneous plates. 

Harris# G.3.. [2 1 studied the b^uckling of xaniaxially 
and biaxially loaded orthotropic plates using Ritz variational 
technique. 

Warren and Norris f 3 3 that for the laminates 

containing laminae of identical orthotropic properties and 
thicknesses. It is possible to orient the layers in such a 
way that the resultant elastic stretching behaviour is 
isotropic. According to them# to achieve this result# the 
total number of layers must be^ three or more and a typical 
Kth layer of a N~iayered laminate must be oriented at an 
angle = {iC-l)/n« However this laminate is not isotropic 

with regard to bending- stretching coupling# 

Barthalmer [_ 4*} suggested that for uncoupled# orthotropi 
laminates made of laminae of identical thicknesses and material 
properties# it. is possible to achieve inplane orthotropy as 
well orthotropy in bending by selection of suitably ply stackir 
order. The fibres lie synraetrically about the orthotropic 
axes of the laminate. 
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Aiaoartsumyan Fs 1 recognised the effect of coupling 
bet*K?esn inplane stretching and plate bending* This coupling 
effect is the major difference between elastic response of 
an arbitrarily laminated end homogeneous plates* 

Hois ten 3 analysed the buckling response of ortho- 
tropic plates with three simply supported and the remaining 
edge free by developing an approximate equation for the 
biickling load. 

Sharma# Iyengar and Murthy J have shown that it is 
possible to eliminate the coupling effects in four and six 
layered antisyiwaetric cross ply laminates by a suitable choice 
of thickness ratios of laminae* Thickness ratios are shown to 

be independent of material properties* 

% 

Jones 1 analysed the bxjckling and dynamic response 
of antisymmetrically laminated angle ply rectangular plates 
using closed form approach. 

I*(usner and Stein [9 ] used a finite difference energy ■ 
analysis method to make parametric studies for angle ply 
simply supported and clamped graphite - epoxy plates having ; 
a large number of alternating plies {+ 9 )* 

Paul# Jensen# PI inch examined the effects of 

inherent mechanical ..couplings exhibited in 'f tilly anisotropic [ 
(i.e* unsymmetrie) graphite/epoxy' laminates on the buckling 
loads and mode shapes* The results indicate that the coupling 
especially those which relate stretching and bending behaviou 



ca\ise out of plane deflections prior to biKkling and 
reduce the btackling load significantly. The analytical 
results also show that the mode shapes exhibit twisting 
due to mechanical coupling. 

Klaus Rohcver fllj has analysed the plate with 
different botmdary conditions edge ratios and fibre angles 
under tiniaxial as well as biaxial normal stresses to see 
the bending-twisting coupling on the buckling load of symme- 
trically stacked plates. An indicator is established which 
allows to approximately determine the btckling load reduction 
with the aid of stiffness coefficients alone. For plates 
mder compression it allows to decide prior to the buckling 
analysis whether or not it was necessary to include coupling. 

Whitney and Leissa |*i2j have given closed form solution 
for simply supported antisymmetric cross-ply plates having 
S-2 type boundary conditions and simply supported antisymmetric 
angle-ply laminates having S-3 type boxandary conditions. 

Jones jjl 3 I discussed the symmetrically laminated odd 
regular angle ply alternating in the sequence ’+©, -6, +©*»*,+©, 
Here the bending -twisting coefficients and 

largest for small number of plies and. become smaller relative 
to other coefficients as N is increased* ; 

wittrick [,14j achieved exceedingly complex solution in 
symmetric angle plates for infinite striiDS in SSSS# SCSC, SPSF, 
CCCC boundary conditions* i 



o 

Reissner and Stav-:,'"-:; p-'U included bending-stretchlng 
effect in entisyiraueti'icel angle ply plates. He showed that 
the buckling load is r-i-iuced due to coupling which is the 
strongest when cr.l_, c. aviall number of plies are used and 
decreases as N Im::- ~:xsesm 

Ashton jlcl suggested an approximate theory involving 
unsymraetric 1 However it is found that this theory 
was not accur.t - fcr all problems. 

Barwcy p. 7l has made an attempt to achieve specially 
orthotropic behaviour of antisymmetric cross and angle ply hybrJ 
composite Ismincte-s taking fibre orientations and thicknesses 
as design variables. He has shown that at optimum point lamina'! 
tends to decouple. He used Galerkin's procedure for solving th? 
b uckl in g pro bl era • 

Sharma and Iyengar [_13j developed an exact solution 
procedure incorporating the coupling effect by combining the 
three equations of classical plate theory to a single 

eighth order partial differential equation. He took levy's 
type of boundary conditions and converted differential : 

equation to eigh'fct order ordinary differ^itial equation and ; 
then to a eighth order algebric equation* He obtained eight i 
by eight determinant after substituting the assumed general 
solution in eight homogeneous boundary conditions along the 
rectangular edges. The uniaxial btackling load has been found ; 
by a trial and error procedure after equating the whole qf 

determinant to zero* ( 
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Whitney ' 1 , tiainc; the? Foui^iov s-:.i'ie 3 approach studied 
the effect of b :y conditions upon tlio frequency of free 
vibration of an unsyrcnetri'? laminated rectangular plate* 

Bert end Mayberry rauc.ied the free vi3o.ratlc>n of 

unsyirraetric ally iarainct. .. „ inotropic plates with all clamped 
edges. He used Raylc: ; '-its procedure for finding the solution 

r *1 

Lin and King ; "'ij studied the free- transverse vibrations 
of tinsymmetrically Ic .’.L'tsd plates. He solved for various 
boundary conditions using Bolotin asymptotic method. The 
method is applicable only if the governing equation can be 
satisfied by products of trignometric functions in spatial 
coordinates paralleling the edges of rectangular region. The 
approximation to an eigen function is expressed as the sum 


of a generating solution and several corrective solution 
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Reddy S. Chao 24 ^ made a comparison study of closed 
fox® and finite eleasent solutions of thick antisymmetric 
angle ply laminated plates. He solved the buckling problems 
with all edges simply supported and all clamped boundary 
ccnditions. ■ ■ 

t .3. LIMITATICMS OF EARLIER METHODS 

(a) Galerkin*s technique by Saxena I 23 ^ requires consi- 
derable computational time. The programming part is 
length and chances of committing errors are also more. 

(b) Eayleigh-Ritz technique by Kamal & Duxvasula 223 
has the same limitations of 6alerkin*s technique. 

This technique is also not that accurate when compared 
to Galerkin*s technique as the displacement functions 
to be selected satisfy only the Geometric boundary 
conditions. For better convergence a large number 

of series terms are required, 

(c) Finite difference method by Hoingner « Stein jl 9 } 

requires the discretization of the plate into various 
mesh sizes. The accuracy of the solution depends on 
the mesh size and smaller mesh size leads to large 
number of algebraic equations solving, 

(d) Finite element method by Reddy & Chao [ 24 3 has less 
limitation but it cannot be efficiently applied to the 
analysis as it requires lot of computer storage capa- 
city depending on the number of elements used. 



o 

o 


Bolotinas asymtotic method employed by Lin and 

!». -*1 

King -as considerable computational difficulty and 

the ...letliOu is suspactable when one or any of the rectangu- 
lar sages of the plats happened to be f re ' » 

f. Closed form solution technique b'y To:-esC>3j Hwted to 
simply 1 ; ^niitions of comiDOsite ia:r.;'.ae::. ss. 

g» Exact sclution tsclinique developed by Sharraa fioj 
requires the solving of by ..light determinant by., trial 

and eri"Oi' procedui'e in thin; ;'...lsyimiistx'iG angle ply laminates. 
In optirriisatioxn probie-a the analysis has to be repeated for 
parameters like fibre orientations# volume fractions etc. 
Overall coraputer- time will be increased at the final stage. 

There are many more methods available in the liter- 
ature for the analysis of composite plates but 'hardly any 
author has used s-uccessive iteration metiiod for the analysis 
of composite plates. Keeping in view oveK> tte limitations of 
the conventional methods# the application of successive 
iteration method has been investigated in detail in the 
present study. 

1.4 , Scope OF the presemt work* 

h study has been made using sixscessive iteration 
technique to the analysis of thin decoupled and coupled ' 
antisymmetric angle ply laminates. The method can be 
applied to the cases to whicih bevy’s method is applicable# 
i.e, boiandary conditions and loadings are such that# 
governing set of partial differential equations reduce to a 
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set of ordinary differential equations. I'his puts some 
restriction on boundary conditions and nature of the 
larainato* Tao of the edges which are opposite to each 
other must be simply supported and others can be of any 
combination. Levy's type of solution is possible for 
antisymmetric angle ply laminates* Antisymmetric laminates 
are considered in the present investigation as they posses 
greater flexibility v^ith respect to d'" '-ign parameters* 

In this study (^. * *©{-0 j© j”©, .]) antisymmetx'ic 
angle ply lairsiaate is chosen as optimizati-.on studies [^173 
reveal that (••• -e,| ©2 | 1 -©2 1 • • • ) anti- 

syrranetric laminate does not lead to any significant increase 
in load carrying capacity* 

The method presented can be applied, in principle 
to all types of plats problems «• bending, buckling and 
vibration fox* which the number of independent variables in 
the govexiiing equations can be redixied to one* Out of these# 
relatively more difficult to solve are eigen value problems* 
Detailed investigation here. is restricted to eigen value ; 

problems* Biaxial buckling and free vibration problems have ; 
been investigated in detail. .The necessity to undertake 
such a study in thin 'coupled composite laminates is to p 

extend the application of successive iteration technique to ! 
the analysis of higher order shear deformable decoupled compos: 
plates if the teehniqi:^ works well* ^ j 
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A coinposite plete has its length to thickness 
ratio (a/h) more when corapared to metallic plates. In 
the plate analysis the; effect of shear deformation is felt 
when (a/h) is less than 20.- For a/h c 20, the effect of 
shear deformation cannot be ignored- Higher order theory 
developed from variational principle by J.N» Reddy [25 j 
leads to a set of 5 sets of governing partial differential 
ecfcations in five displacement variables u, v, fx and^y. 

In antisyriimetric angle ply laminates it is possible to eliminate 
the decoupling coefficients g, ^26 suitably 

manipulating thickness and stacking sequence. The problem 
thus reduces to the solution of a set of 3 governing equations 
v;hose basic nature is the same as that for coupled laminate 
using classical plate theory* 


In Chapter II of the present work, expressions are 
given to generate laminae properties from fibre and matrix 
properties. The governing differential equations for buckling 
and free vibration of rectangular laminates plates are presen- 
ted. The boundary conditions are considered with x edges 
simply supported and y edges simply supported (SS), simply 
supported-cl amped (SC), clamped-clemped (CC), cl amp ed- free (CF), 
simply supported-f ree (SF) and free-free (PF). The development 
of successive iteration rnethod to the analysis of composite 
plates has been discussed. Analytical expressions are given, 
which ensure decoupling in six layered antisymmetric angle 
ply- laminates. 



In Chapter II 
biaxial buckling and 
laTiinates, using this 


li, an attempt has been made to scZ.ve 
free vibration problerxis of decoupled 
rnethod. All the boundary conditions 


given above are considered* 


In Chapter IV« investigation has been made to appl]; 
the method to the buckling analysis of coupled laminate vvit 
edges simply supported and y edges clamped condition. Var; 
computational model schemes have been tried v;ith three and 
dependent variables* 


h X 
oils 
two 



CHAPTER -• II 


GENERAL ANALYSIS OF COMPOSITE PLATES 

A s'cudy has been mads to apply successive iteration 
technique to the analysis of composits laminates* As 
decoupled antisyir^Tie trie angle ply laminates are involved 
in this study# expressions for decoupling thickness ratios 
are presented in this chapter* 

2.1 ASSUMPTIONS IN THIS COMPOSITE PLATE THEORY 

The laminated plate analysis is based on the follow- 
ing assumptions as sifted ty 

a* Each layer is orthotropic# linearly elastic and of 

constant thickness* 

b. Thickness of the plate is very small compared to 
its length and width* 

c. Body forces are absent* 

d. The displacements u# v# w in x# y and z directions 
respectively are small compar-ed to the plate 
thicknesses. 

e. Strains g’.,# and are small compared to unity* 

.A i 

f* The plate is in a state of plane stress* 

g* Transverse shear deformation -is neglected* 

Rotory inertial terms are neglected. 
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2.2 ANALYSIS 


Consider a N-layered laminate as shown in Fig. 2.1. 
Each layer, made of continuoiis fibre reinforced material# 
considered to be orthotropic with uniform material proper- 
ties ^ 2 .- 32_2 and. Yi2* stibscripts 1 and 2 refer 

to the axes along and perpendicular to the fibre direction 
respectively as shown in Fig. 2.2* For such a laminate 
stre'ss-s train relationship of iayev* IS expressed 

by Jciftsslisj as 



— 



0 



0 


^2 



^12 



^ J 


C 2 . 1 ^ 


where are the reduced stiffnesses 
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^“21 ®1 


(1 -Yi2r2i) 
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■ ** %2 


®22 


# *■ ^12 ^21^ 




The Foisaon*s ratio can be obtained from the relation 





If -the fibres in a lamina lie at an angle Q to the x-axis of 
reference plane. Fig- 2*3, then the transformed stiffness 
matrix for the K-tli lamina is given by. 
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+ Q 22 - Sin^© Cos^© + Qi 2 ^Sin^« t Cos^©) 

^22 
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4 2 2 4 . ^ 

Sin © + 2CQ|2 ♦ 2Qgg) Sin © Cos © + Q 22 Cos © 
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- Q |2 * 2Qgg) Sin© Cos^© + ^^12 "** 

°26 


^%i 

- 05^2 ~ 2Qgg)Sin^e Cos© + ^Qi2~^22 ^ 

^66 

m 

CQli 

‘‘‘ ^22 “ %2 * 2Qgg)Sin^©C03^© + Qgg(sin^© + Cos 


(2*5) 


Frc«B eqns* (2*5), .it is evident that in lamina with 
fibre oriented at to x axis having identical properties, j 
all the elements remains the same except and Q 2 g wh4 

differ only in sign* 

Mat haaa tic al 1 y, 

** ** i*l,2 J 



15 


Using KirchQff's laws thin plate hypothesis with usual 
notations# the following constitutive relations 
obtained W 
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(2.7) 


vdiere# 

a ‘*'*^*yy \y ~ ” ^^'xy strains and 

N^r Uy# Njjy# My» are the stress resultants. 

The exterosional coupling and bending stiffnesses which 
are expressed as f^n* and D| ^ respectiTely are 

^ij 

■' 

®ij 

hi 




(2.3a) 


( 4 - 4 . 1 ) 

(2.8b) 

^ ^ij 

<4 * *Li> 

(2.8|f)C 



For the antisyniraetric angle ply laminates that are used in 
the present sttidy 

*16 “ ^6 - ®16 “ ®26 “ %6 - “26 “ 0 

using strain displacement relations# eqn* (2*7) can be rewritten; 
for antisymmetric angle ply laminates as 

** ^1 * ^12 ^ %2 '*''yy 

®y '*■ ^^22 ^'y " %2 ■“ ®22 '^'yy 

N « Agg (u#y + ir#jj) - 2Bgg w# 

• (2*®) 

^ * ®11 ^x ®12 "^'y * ®ll ^*xx *" %2 ^'yy 

^y “ ®12 ^X ®22 "^'y * ®12 ^'xx ■* °22 ^'yy 

^y * %6 ^^'y ■*“ “^^x^ * ^*Xf 

Puting £ ** x/af h •* y/b and defining length to total thick- 
ness ratio as r « a/h# aspect ratio as p « a/b eqn. (2.9) 

becomes ' : ; 

“-£ = < r > "•£ + ^2 ‘ e > ’^' 1 , - ® ll '^> ’'• m " "'rin 

“*1 “ ^2 “■£ ^2 '''►1 ■ *'' 6 £~ ® 22 < l >^ ’'' H'l 

»£>, ’ ^8 ”r) '*"1 + <r> - ^®86 
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« = h (i) u,,; + Bj 2 ’"'££■■ '"'ll* 

^ IT 

M = h (ii 2 (Ij u,^ + (§) v,^ - 522(S> '"'fet- ° 22 ‘?>^ "'nn) 

♦ 

M s h (Sgg { (|) u, »^ + (i) v,£ ) • 2Dgg i^) ) 

r 

(2*10) 

2,3 GOVERNING EQUATIONS OF ANTISYfffiETRIC ANGI 4 E PEY EAMINATSS 
Governing equations for biaxial buckling and vibration 
as Statadi by 3b»nes £*2Q as 

*11 ^xx *66 ^yy ■*■ ^* 12 '*'* 66 ^^'xy * ^%6 '^'xxy‘*®26’^'yyy * ® 

( 2 , 11 a) 

^*12‘*'*66^ u,^y + Agg '^*xx * *22 ^'yy * \&^*xxx * ^^26 '^'xyy * ® 

( 2 . 11 b) 

%1 ^Jxxxx ^^°12 ^^66^ '^'xxyy ^22 ^*yyyy “ ®16^^^xxy”*'’'^'xxx^ 

-»26^^yyy + ^’^'xyy^ ^"'xx '*'yy ■^^'^'tt “ ® 

( 2 , 11 c) 

Figure 2.'4 sHovv© the loadiifig and geotm+py oF the laminate • 
for free vibration probl«a# in the above equation, set N^ *» Ny « ( 

and for buckling problems setPw, •* 0, Eqn, 2,11 can be 1 

rewritten in non-dimensional form as i 

1 

* 66 P^'"%^ <^ 2 +^ 6 >P'^ "'eeri- “ i 

( 2 , 12 a) I 

I 

»^+ *22 v^,,+ ^66 •''eg- 382 * W'e^rj 
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'3ii6 ^26 ^ ^^enn" ■“ 


*** ^^^12 3Dgg)p ^^*c.c:nn^ ^2 7 P W<iTi'>»>l'*' ^^ * '££ "^ ^ P ” ^’^*00 ^ "* X (w) 


'66 


'8£r>r^" "'Sa 


“*nn 

(2,12c) 


wJ-r-;: .re# 
X ss 


li„ b" 

WM^WMWMWMIt 

Ej h^ 


and X 


p b^ (30^ 
''' ■■''TY"”” 

^2 ^ 


N 

K » 


H. 


X 


For biidcling problem “ 0 and for free vibration problem A » 0 
in (2 ,120). 


2,4 BOUNDAKJf CONDITIOHS OP AHUISYflMETRIC ANQIiK PLY LAMINATES 
The governing eqns, (2.12) are solved for the follow- 
ing botmdary conditions as -furnished by Tones Dal* Boundary condition 
on X edges in all the eases under investigation are sin®>ly 
supported and on y edges six combinations are investigated 
on X edges# 

SIMPLY SUPPORTED-S IMPLY SUPPOIYB© 
at E = 0# 1 

m = w * M^ = » 0 

2,4*1 SIMPLE SUPPORTBS-SIMPLI SUPPORTED (S-S) 

at n « 0»i (2*13a) 

V ®* W s» ii, = M sx Q 

w tr\ 




2.4.2 SIMPLY SUPPORfED-C.<AMPED Cs~C) 


at 


0 






V " 

= W = 

II 


= 0 


(2i 

»1 3t) 

at 

‘1 = 

1 





u 

« 

as w 

8 

II 

0 



CI^PED-CLAMPED 

(c~c) 




at 

m 

0,1 






u 

ir 

= w 

8 

w# n 

0 

(2, 

.13e) 


2.4.4 CLAMPED-«PREE (C~P) 

at ^ = 0 

ti = V =* w *» * 0 (2«l3d) 

at n = 1 

M = M =. 2 + p = 0 

2.4.5 SIMPLY SUPPCaiTES-PREE (S-F) 

at s 0 

V a w ■■= « 0 

at »^ =» 1 

“n “ % = ”n “ ^ “n -n “ “ 


^0 

2,4.6 PREE-PREi: (P-F) 
at *1 * 0#1 

{2,13f) 

2-5 LEVY'S SOLUFIOK PROCEDURE 

The boandaxy conditions considered in the present study 
ars simply supported on x edges. The displacement, functions 
u, V, w can be assumed in a form which when sxibstitnted in 
eqns. ( 2 . 12 ), reduce these equations to a set of ordinary 
differential equations. 

Assuming, 


^ m 

u(n) 

sineatfe where m :nn 

and m is a natural nuwto^ 




C 2,14a) 


¥ (n) 

cos^Q. 

(2,14b) 

w ^ 

w(*^) 

sincig 

(2 ,140) 

Siibstitiitiiig 

these 

expressions in (2.12) 

one gets the following 


set of ordinary differential equations 


f I I i lit I 


e, w"^^ e« w 

il*. .m 

+ fcj w 

+ e, ¥ 

4 

t OgV 

+ ©g u 

•i-e^u 

(2,15a) 

. m ■ 

■ It 


» 


(2,15b) 

©g w t «9 w + 

®io ^ d 


®12 ^ 

* 0 

f 1 # 


1 

1 



®13 ^ ®14 

W 4* ^ j g 

+ “ts 

u, + 

^17® "* 

0 


(a^isc) 
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wher® 


m 

°22 ^ 

4 

P ; 

£> sa 

^2 

1 -(2 

(I>3^2 2Dgg)) 

2„2 

P 

+ Kp'^'A ? 





' %1 

^ p2 ^ 


,, 2 . 

ir # 

i 

m 

3a26 

oC IT 

p^# 

®5 

* - ®6 

WSi 


mj » 


. 2 
of ** 

rpj 

®8 

* -3i2gC<p^? e. 

33# 

9 


®io • 

’ ^2 

2 

p r# e 

'll 

7 

« — o( " r? 

e. ». 
12 



®13 “ *"^26 ®14 * ©is « + Agg)(~o() rpi 

®16 ** ^6 ®17 * ”'^11 ^ ^ (2*15d) 

It may be noted that for free vibration problems 'X=i o and 
for buckling problems Xf = 0 in the above equation* 


2*6 DEVEbOPMEOT OP THi- SObOTIOli TECHNIQUE 


Writing eqns. (2*15) in the form 



/ • * 

« f Cw # 

w* 

i 

1 

# V# 

1 1 1 

1 1 

V 

1 1 

«, g(w # 

w* 

ITt 

I 

u ) 


t 1 

• ,1 

# 

1 

1 


m ; 

1 

# 

W # 

r t 

u) 


(2*15e) 

('a.lSf) 

(2*15g) 


and integrating eqn* {2*15e) four times# eqn* (2*15f) twice* 
eqp* {2*5g) twice one gets expression for w#rjr|# w*#^ # v*»^» 

V* u#*^ # u* These relations will involve eight constants of 
integiration whioh can be determined from the eight bO'Undary 
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con<iitions. Starting from some inttial.l 3 r chosen functions 
u# V and w it is possible to develop an iterative procedure 
to determine the solution* 

ite the governing equations and associated boundary 
conditions are homogeneous# displacements are indeterminate 
with a scalar multiplier^ on a can arbitrarily fi^ a normali- 
sation condition say w(i/2) « 1* This condition deter- 

mines the eigen value \ at each stage of iteration* However# 
when the solution was attempted for laminated plates iteration 
does not converge* It appeared that numerical instability 
occurred because of the eigen value assuming negative value* 

A detailed investigation# starting from simpler problems# was# 
hence# carried out* 

The proced\ire was investigated for column-buclcling 
problem with both ends. 1. simply supported and 2* clawqped 
conditions* A normalisation condition of w(l/2) ** 1 was 

used for finding the eigen value ^ . Convergence was achieved 
in the iteration and the eigen value matched closely with 
the e3cact result* 

The procedure was then investigated for isotropic 
plate buckling probleu where lievy*s solution procedure is 
applicable* Problems with x edges simply supported and 
y edges in 1* simply supported, 2. clamped and 3* clamped- 
free boundary conditions were been studied# using normali- 
sation condition w{l/2) * 1 * Convergence of iteration 

was observed In the first ease of boundary conditicm for all 



c3 

aspect ratios between 0*2 to 4.5* But in tbe second and 
third case of boundary conditions oscillatory tendency in 
convergence with eigen value becoming negative in the 
iteration process was obaej'vedjfor many aspect ratios* 

It appears that X mxast be positive at each stage of 
iteration to converge. To ensure this# a separate expre- 
ssion for X is necessary. I'his can be done by obtaining 

I 

expression in the form of a Rayleigh .s quotient by using 
energy principles* 

In the mame isotropic plate problem ‘cases investi- 
gations are carried out with Rayleigh's quotient X Instead 
of the normalisation condition* The oscillatory tendency 
was not noticed and the convergence of the iteration was 
achieved# agreeing well with the exact result for all the 
aspect ratios. 

For laminated plates, the need of Rayleigh's quotient 
expression is thus confirmed. 

2.7 KXPRESSICM FOR THE SIOM 

The total potential energy expression for the problem 
under investigation is ofctalned whHney D‘>-3 

. o o . ** 

+ (u#y + ** ^*X3K 

t 2u#^ ^'yy '*^'yy ^y#'‘^*xy^ 
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+ D22 * “66 ~ " *'y(«'y>^ 

«pco'- (w)^)3 dx dy (2.16) 


For buckling prob.lerri 0 and for free vibration 

problesn » N «= 0» using the variables and , 

Y 

Eqn. (2.15) crn be rewritten as 


1 1 


T ' ,v27 - ...2 — 22»» 

r + 2 Aj^ 2 ^22 ^ ^ (v^.q ) 


o o 


■*■ He + 2pxi,^ v.£ 4 (v,g)^r^ ^ 


^ 3 i 

+ v,g_ w,g^^+ 2u,^ w,g,j^ pr ) - 2 B 2 g (u,r^ rp + v^ \ 

rp^ + 2v,i^ w,£^p^r) + (w,£^)^+ 2d^2 

+ ^22 + 4Dgg p2 - )\((w)^ + Kp^(w,fj)^) I 

( P^ (w)2) ] d£ d»^ (2,17) ’ 

It may be rementioned that H " ^ buckling problems and 
A =« 0 for free vibration problems. Using expressions (2.14) 
for u# V and w in eqn* (2*17) one gets 

p * T -T ^ (uv ) -1- A 22 Cv )^ 

+ Agg (p^ (u - 2 podvLV) + v^«<^)r^ + 2i|g((m w) proc^ 

• Cw)c<^ r - 2(uw })o("^ pr + 2 B 2 g (~(u w ) rp" 





**2 *’2 ••4 2 «“'22 

+ Cv w ) rp «: • 2 ( V w )<^ p r + 5. , a< (w) • 2D^« p oC 


(ww**) + D22 (w‘*)^ + (w’)^ + ^[-P^ ( ^(w)^ 


66 


+ (pV)] d£ dq 


( 2 . 18 ) 


7\ 


C^il ■*■ Pi] " ' 


( 2 . 19 ) 


where C’^jl the internal strain energy component of the 

plate. [^ll r)Otenti;il. energy component due to 

the inplane load and Xf [Q 2 I kinetic energy component 

due to vibration. 

For a neutral equilibrixan total potential energy ig 
a constant. This implies that can be obtained by setting 


Xr * 


■P 


Doing this X « " / pi} 


(2 .20a) 


Natxiral frequency is determined from the fact that 
the maximum potential energy and maximum kinetic energy during 
an oscillation are equal. This implies 


% ” Pal 


(2.20b) 


Aa already mentioned the application of this technique 
is studied in coupled and decoxilpled thin antisymmetric angle 
ply laminates* Since decoupled laminates are used in the 
analysis separate analytical expression to ensure desoupling 

is 3r®gttir#do 
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2.9 EXPRESSIONS FOR DECOOPLIiNS THICKNESS RATIOS 

Consider a si* layered antisymmetric angle ply laminate 
as in Pig* 2,5. In such a laminate for every lamina oriented 
an angle o above the mid plane there exists an identical lamina 
placed at an equal distance below the raid plane but oriented 
at a negative angle O* One needs to make and zero for 
decoupling between inplane displacements u# v and transverse 
displacement w# For given material properties and fibre 
orientations of the laminae it may be possible to malce 
and I 26 zero by suit^ly selecting the thickness ratios 
the layers of a laminate. Proceeding in a similar way as 
Sharme. [,17"] and Saxena £24 ] • Using eqn* (2. 8b) 


0 ) + Q. 


((t2+t3) 


t^) + 



Eqn. (2.21) becomes 

Pij + ^ i" ^2^^ — 1 ) + ( (1 + + R 2 ) •(It ^ 2 ) ) ■» 0 

( 2 . 22 ) 

Let (1 + ^ 2 ^^ “ (l + R^tR 2 ^ * ^2 

Plj * '^1 - * 

solving e<3[oiCirx') 


0 





./ 


1 + A 


R- 


/XT' ” ^25^ 


6 ””^16^ ^ ■i'A (2*23) 


where A - iPl6 " ^26'' ^'^26 ' °‘l6> 


Plj‘^6 


for R, >0? A> 0 and for Re > 0; ^ i 

®^ 16^26 


For .iOyeneous laiiinatas with alternate + Q aiid -9 ( Pig® i-s") 

^ j_j =a 1 anc, <x^j » -I* Substituting in eqn.Ca-i^-) 


R * _(B + 1) + -'aRl + IRx 


(2,24) 


Any coirtsinations of Rj_ and R^ satisfy the above 

relation, decouple the laminate. 

The analysis is much simpler in decoupled laminates, 
as the number of governing equations are x-educed to one with 
this ratios. The application of the successive iteration 
technique can be studied readily in decoupled Ismiinates. 
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2 UN IDIRECTlUNALLY REINFORCED LAMINA 












CimPTgR ■ III 


ANALYSIS OF PECOUPLED LAMINATES 

In this chapter an attempt has been made to analyse 
decoupled six layered antisymmetric angle ply laminates using 
successive iteration technique* Biaxial buckling problems 
are studied with x edges simply supported and y edges in 
combinations of (l) simply supported# (2) simply supported - 
clamped Cs-C)» (3) clamped (C~C)# (4) clamped free (C-F)# 

(5) simply supported - free (s-P)# (6) free (F-F) boundary 
conditions* Free vibrations problems are studied with x 
edges simply supported and y edges in combinations of 
(1) simply supported (S-S)# (2) clamped (C-C)# (3) free (F-F) 
boundary conditions. 

3.1 BIAXIAL BBCKLIWS FROBLlflS 

3.1.1 GOVERNING EQUATION 

Since SjLd '®26 R*ade zero# the eqns. (2.8a#b#c) 
get decoupled. The eqn. (2.86)# then reduces to an equation 
with a single displacement variable w as 

w'*’ 4- ®2 *** ##0 (3.1) 

where 

• ^22^^^ ®2 “ + Kp^ X # o<^ 
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Eqa, C3,l) can be rewritten as 


w 


IV 


Cj^w + CgW* 


(3,2) 


where 


w I X f a. 


D 


1 . , 2 

” (~) f z. 

Jb^ 

22 


•z— ; 
“22 

2(Dj^ 2 + 25gg) 


*# 


K 


D 


22 


(3,2a) 


W** ' a* 


Integrating the eqpi. (3,2) successively 

h h 

Cl J i/di^ 4^ Cp I rj 

o « 


(3,3a) 


w“ » 


W ™ ffla f*** 

*** \#> gj 


w « c- 


n 

8 r 


.^-2 


I' i -4“ C 2 W ^ '^i^) + K- 


o c 


!) ^ 


(3,3b) 


J 1 J wdrj^ -i- Cp jwdr^ + *3.^™ + K2r| + (3,3c) 


0 o a 

Vj 7 /' 

0000 


n 3 

wd + C2 ^ ^ wdl^^ -t 






(3.3d) 


where K^# K 2 » Kg and K^ are the constants of integration to 
be determined fton the boundary conditions. 



3.1.2 RAYLEIGH’S QUOTIENT 

Governing eqns. (2*i2a»b,c) get decoupled for zero 
values of and B 2 g. Eqn* (2.12c) can be repeatedly- 
integrated to determine w and A , Eqns. (2.12a) and (2.12b) 
are the coupled equations ds'carraining u and v* However these 
equations are homogeneous and the associated boundary condi- 
tions are also homogeneous. In the absence of any eigen para- 
meter the solution of these equations is u =* 0 and v « 0. Using 
these zero values for u arid v in eqn. (2.18). One gets the 
f rillowiu,^ expression for \ from eqn. (2. 20a). 



3.1.3 (34S) BOmiDARY CONDITION CASS 

The bo\mdary conditions are obtained after substituting 
the displacement relation (2*l4c) in eqn. (2,13a) as 

w 4i w* » 0 at « 0itl {3.5) 


Let the definite integrals be defined as 
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Using eqn. (3*5) in eqns. (3*3. b) and <3«3d) the constants 

of integration are obtained as# 

Kj^ w - (C3_<l2 t C2 w{l)>; iCj *s 0; Kg « - g - ^4 + ^^4*® 

(3*7) 

3*1.4 STEPS IK APPLZIMS SUCCESSIVE ITERATIOK TECHNIQUE 
The main steps to apply the siiccessive iteration 

teclmiqne are as follows i 

1} For the gir^en fibre and matrix properties Ej^, £ 2 # 

Gj _2 ^12 a te Cq 'J and [.dlnafcrices \is ing 

eqn* (2.4) and eqn. (2.5) respectively, 

il) Malcii^ iise of the decoupling thickness ratios Rj and 
R 2 calculate the percentage thicknesses of the 
1 ayers • 

iii) For the specified geometry and thickness# compute the 
coefficients of the stiffness matrix 

iv) For the given aspect ratio Cp) and mode ntcnber (m) 

calculate the constants ^ 4 * 

v) Select some initial function w(if)) and an initial 
value for X • Choosen function need not (say wCiq) « 1) 
satisfy any bomdary conditions, 

vi) Find the constants c^ and © 2 * V^ven boundary 

conditions of the problem evaluate the integration 
constants K*# JC># K, and K*, For (s-S) case find 
the cons tana ts from eqn* -'( 3 # 7) ♦ 
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vll) 


viii) 


ix) 

k) 


xi) 


Evaluate the functions '^n+l' '^n+1 eqns. (3 #3) 

using where suffix n denotes the iteration n\amber. 

normalise '^n+l ^n+1 raaximtm value 

of Wij. 

Find the eigen value ^ f rom eqn* (3*4) ♦ 

Transfer the control to step (vi) and repeat the 
process until | ( \ ^^ . X_^)/ ^ g, , a small 

nuiTiber (say £ 0.001). 

For various values of mode nuoibers# find the eigen 

values ^ and the minimm ■. out these is the 
buckling load. 


3.1#S CONSTANTS OF INTEGRATION 

For the remaining various boundary conditions at 
the y edges# the constants of integration are found in 
a similar way as in SS case. I»et the constants be defined 
as# 


a. 


D. 


12 


"22 


«3C. 





^Ii2-L±hs}, 

°22 


®3 “ °1 ^®2 * *1^4^ * °2 ('*<!) + ®ji “ °l(% + ®2‘*3^ 

t C 2 (w* (1) + ag «* Cl + *» (1 + — ^) I 

a,^ a (1 + ag W 4 * (3*8) 

then for various boundary conditions at y edges. 
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SIMPLY SUPPORTED-eLAMPED CS~C) 

» {3 (q^ - q^) + C2 (q2 *“ 'll))? ^ 0 

Kg ** (0*5 qg *• 1*5 ®2 ^ 0* 5 q^^ •“ I •) 5 qg J t K^ ** 0 

(3,9) 


CLAMPED-CLAMPED (C-G) 

^^2q^ - 6q^) + Cgdaqg * {2qg-6q^)4-C2 CSq^^-Sq,) 

Kg « Of K4 » 0 (3.9b) 


CLAMPED-PREE (CsF) 




(ag -I- Bj Kg) ^ ^ ^ ( a^ ag - aga^) 

^5 (aga-^ - ajag) 


•f Kg**0f K^^O 


{3.9<j) 


SIMPLY SUPPOKTED-FREE (S-F) 

(ag + ®J^Kg) 


WT 


a. 


I Kg • Of Kg 


Ca^ag 




(a^a. 




I K4 «* 0 


(3.9d) 


FREE- FREE 

** -agiCgf Kg 


2 

2 (2a^ag + 

**■ 


I Kg « 


( 2 ag ag Kg ) 

2aQ 


-2. 

h 


(3.®e) 



3.2 FREE VIBRATION PROBLEMS 


The eqns. (3.2) and (3*3) are the same with the differ- 
ence in value of Cj_ 


a* a- - X,,'? z, 


D 


D 


11 # V 4 

- (■— )S ^6 

22 


O. 


22 




^2 “ 


■rtw* I i 

P 


D 


(3.10) 


22 


Expression for Rayleigh's quotient is obtained in the same 
maimer as explained in section 3.1.2# 


X _ * 

f 


o 2 "w" 4D22P^ (w”)^+ 4Dgg dr| 


..2 


o 


w 




(3.11) 


The bomdary conditions and the values of the integration 
constants are the same as that for the budding problem. 

The analysis has been m^e in (S-S)# (C«-C) and (F-F) 
boundary condition cases with staccessive iteration technique. 

The application of this technique is very similar to the 
buckling problem. 
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T ^2 * 0.25 and CPRP are * 0*667* 

” 0.20. The decoupling thiclmess ratios of ** 

and ** t 2 /'t ^2 *" ^*0 are used in the analysis. Hie analysi 

was carried out in DBC-1090 system at.I.I.T. Kanpur. 

Bo’indairj conditions SS, SC, CC, CP, SP aM PP haive been 
considered in biaxial problems to determine the buckling load; 

A ^ using successive iteration technique, Conditior 

E,j 

CC, SS and FP have been considered in the free, vibration probleu 
to determine the non-dimensional natural frequency * ab —— _ 

Sticcessive iteration technique was studied in detail and the 
eigen values' obtained In the iteration prc&e^s . are presented in 
tabular foiro for differoat starting functions, and mesh sizes for 
quadrature using trapezoidal rule. A parametric study has been 
made to investigate the variation of budding load and natural 
frequency ( with different fibre orientations (9) , aspect 
ratios Car) and buckling load ratios (KK 

3.4 RESULTS AND DISCUSSION 

SiKicessive iteration technique gives very good results, 
agreeing well with the results of other techniques. The present 
technique gives slightly higher values of buckling loads and 
natural frequencie*® compared to other exact solution techniques. 

The reason Is due to the use of energy expression and discrete 
iisplacemaat functions in the analysis. 

Table 3*1 shows the eigen values (X) obtained in the 
iteration process for various aspect ratios (AR) with fibre 
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orientation 0 *= 30^ and buckling load ratio K « 0 m 

1 r.vi.u'nate. A mesh size of 0.01 is used v/ith starting function 
w 1 = 1 and convergence criteria ^ = 0*01 • 

for AR 0.5# 1.0# V.6 and 2.0# the buckling loads are obtained 
as 130*9# 79,36# 87*75 and 79*16* ■ In most of the cases the 
iterations converged in not more than six iterations* '1?he 
solution is obtained quickly and constimes less computational 
time* 

arable 3*2 shows the effect of starting functions# raesh 

O 

size and convergence criteria in GFRP , laminate with 9 ** 30 # 

AR => 1 and K a 0* 

i) 

ii) 

jik^Lii 5 

ir) 


For a starting ftanction w()q) =1, mesh interval = 0*01 
and convergence criteria £ = 0.01# the buckling load 

( Xjj ) is obtained as 87,75. 

For a starting function w^) = Sinnir^and with same mesh 

size of the above case and the buckling load C ) 

' °2 

is obtained as 87*75* 

For a starting function w(q) « 1 and z- 

mesh size *» 0.001 the buckling load (X^ ) is obtained 

as 37.75. 

For a starting function wC>^) =s 1# mesh size 0*01 and 
£ « O.OOl the btwkling load ) is obtained as 
37*71* 
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3.4.1 BUCKIiINa PROBLEfIS 

All the laminates conaidere«a in the present study are 
OPRP except in Pig. 3-7 where CPRP is taken. 

Figure 3.1 shows the \ vs. © curves in CC# SC# SS^ CP* 

SP and FP conditions with K = 0 and AR a 1. !Phe curves in GC 
condition agrees closely with sharma tl83 and in SS condition 
agrees with Jones [l3 3 • As © increases, "X increases# reaches ■ 
a maximum value of 102.3 at © » 43 in CC condition# 88.2 at 
© * 44 in S3 condition and 68.1 at © » 45 in SS condition* In 
all free edge cases# as © increases# X decreases- \ reaches a 
raaxiratam value of 34.2# 32.3 and 30.8 in CP# SP and PP conditions 
respectively at © = 0. It is Observed that|iX(CC) > X (s3),>'X(ss); 
X(CP) >X(sP) > X (FP) in the curves for all values of ©. 

Figure 3.2 shows the vs. © curves in CC# SC and SS 
conditions with K =» 1 and AE * 2 • "Phe same trend as in Fig.I.l 
is observed* reaches a maximum value of 42.1 at © ** 66.1 in 
CC condition# 24.4 at © * 72 in SC condition and 13*2 at © • 7i 
in SS condition. 

Figure 3*3 shows the Kvs. AR cuirves in CC# SC# SS# CP# 

SP and PP conditions with K » 0 and © *» 30°. As the AR increases# 
X decreases in all the curves upto AR • l.S# and then beccwe# 
more or less constant. The maximm £o>- ipc# sc# ss# CP# SP 
and PP conditions are 130.9, 129.8# 114.9# S3.8# 51*8 and 12*1 
respeetivoiy# occurring at m « 0*S* Tlwt curve in OC condition 
agrees with Sharraa j^lS] • 
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Figure 3*4 shows the X vs* AR curves in CC» SG» SS# CF, 

SF and FF conditions with K =» 1 and © « 6o°. The same trends 
as in Fig. 3.3 is observed. The maximum X for CC, SCf SS, CF, 

SF and FF conditions are 68.8, 50.9, 42.3, 21,9, 19.8 and 8.9 
respectively, occurring at AR = 0.5. 

Figure 3.5 shows the X vs. K curves in CC, SC, SS, CF, 

SF and FF conditions with © = 30 and AR 1 • As K increases, 
decreases. The maxirauirt Xfor CC, SG, ss, CF, SF and FF condi- 
tions are 77*8, 67.9, 58.9, 28.1 and 8,8 occuurring at K «* Q 
Cmiaxial loading). 

Figure 3.6 shows the X vs. K curves in CC, SC, SS, CF, 

SF and FF conditions with © = 60 and AR » 2? The sme trend is 
observed as in Fig. 3*5. The FF ciirve reaches almost aero value 
K 0.5. The maximum for CC, SC, SS, SF and FF coi^itions 
are for 62.1, 53.9, 44.0, 23.1, 19.90 and 9.1 occurring at K « 0. 

Figure 3.7 shows the X vs. AR curves in FF condition 
with K *K 0 and AR *» 1 for © 30 , 45 and 60 . The curves agree 
closely with Sharma 18 . As AR increases ^ decreases. 

^ (# « 30) "> (#»45°)>' a 60^) is observed. 

• .4.2 FREE VIBRATICai PROBLEMS 

Figure 3.8 shows the.fXf curves in CC, sS and FF 

condition with t • 30 . As AR increases, Jfx. f decreases upto 
AR « 1.2 in CC condition and' 1.5 in SS condition. For AR -pl.a 
in OC condition and 1.5 in SB conditionix f increases. The 

for CCand SS conditions are 33.9 and 32.8 
i?eiq^ctiveX)f « Xn FF condition as jypi increases rx value 
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decreasest continuously from 32 at M =* 0.5. 

Figure 3*9 shows the T'Xf vs» © curves in GC-# SS and 
FF conditions with AR»»l.o. increases as © increases and 

attains a maximum value of 41 at 9 » 90® in CC condition* In 
FF condition, decreasing trend is observed as 0 increases* The 
maximira rX^ occurs at © =» 0. In SS condition, symmetry is 
observed about © « 45® with a maximxim value of ^ as 25.8 and 
the curve agrees with Jones \jl3 J • 

Figure 3 •10 shows f ® curves in GC, 38 and FF 
conditions for AR «» 2*0. Cx ^ increases, as © increases in CC 

-mi, ^ 3 ^ 

condition, reaches a maximum value of 'i = 82*2 at © «* 90 '• 

In SS condition as © increases^ ^ inci'eases reaches a maximum 
of 36 at © * 84.3* In FF condition d«;rea#es, as © increases 
witdi maximum occurring at 9 « o* 

3*5 GFKERAL OBSERVATION 

Based on the above study using successive iteration 
technique certain general ccaaraents can be made. 

1, Successiire iteration technique gives veigr good results 

agreeing well with the results of other techniques 
in both bUBhling and free vibration problems- The 
techniqw gives slightly higher values of buckling 
loads or natural frequency compared t© other exwst 
solution techniques- The reason is due to the use of 
energy expression and discrete displacement functions 
in the mialysis- 



The computational efficiency is quite high and in most 
of the cases the technique converges in Just 4 or 5 
iterations • 

A mesh size of 0«01» interval is sufficient^ for quad- 
rature usii^ trapezoidal rule, to achieve desired 
accuracy. ^ 

Convergence criteria (| | > o.Ol) is 

^ n 

siifficient for the two successive iterations to 
converge. 

In buclcling prohleas with OC, SC and ss conditions as 
fibre © increases# bxachling load X increases, reaches 
raaximtwa for particular © (optimal) and then drops* In' 
plate with CF, SF and FF conditions X value drops 
continuously as © increases and hence it ,is desirable to 
design plate for zero degree fibre orientation to 
achieve optimization* 

As the aspect ratio AR increases > decreases in all 
the boundary conditions* 

As the buBWling load ratio K increases, X decreases. 
Interaction curves for biaxial loading;.® an be drawn, 
which will' be useful,, for design* 



FIBRE ORIENTATION (0) »- 

VARIATION OP NONDI MENS I ONAL BUCKLING LOAD (X) 
WITH FIBRE ORIENTATION (0) . (K = 0.0, AR = 1 



FIG. 3.2 : 


FIBRE ORIENTATION (©) *- 

VARIATION OF NONDI MENS I ONAL BUCKLING LOAD 
WITH FIBRE ORIENTATION . (K * 2.0, AR = 2 




V;iTH ASPECT RATIO (AR) . (K = 0.0, © = 30) 



FIG .3.4 : YARIATIOJ OF NC^IMENSIONAL BUCKLING LOAD (A) 
V.’ITH ASPECT RATIO (ARK (K = 1.0, © a 60.) 



NON DIM3NSICMAL BUCKLING LOAD i» *- NON DIMENSIONAL BUCKLING 



FIG. 3. 5 : VARIATION OF NONDIMENSIONAL BUCKLING LOAD (?») 

WITH BUCKLING LOAD RATIO (K) . (AR = 1.0, 0 = 30°) 



FIG. 3. 6 s VARIATION OP NCM DItffiNSIONAL BUCKLING LOAD (X) WITH 
BUCKLIIK5 LOAD RATIO CK) . (AR » 2.0, « = 60°) 




PIG. 3. 7 : VARIATION OF NON DIMENSIONAL BUCKLING LOAD ‘(A) 
WITH ASPECT RATIO (aR) . (K = 0.0, AR = 1.0)' 



FIG. 3. 8 * VARIATION OF NONDIMENSI ONAL NATURAL FREQUENCY ( 
WITH ASPECT RATIO (AR). (8 * 30°) 




7able 3«1 t of Aspect rsfetos ia Decoupled Laminar^e 





10 

m 

m 

m 






m 

CM 

CM 






« 

# 

« 

• 

« 




« 

in 

tn 

OD 

H 

o 




S 

o 


CH 

O 

o 








H 

H 



o 


o 

o 

m 





• 

w 


m 

T*t 





CSI 

i 

9 

# 

« 





J 

e 

m 

m 

01 







CD 

n-- 

IN' 
















o 

o 

O 






irt 

« 

• 

o 






1 

w4 

tHi 







s 

CM 









iHi 

ri. 

fM 







CD 

r4 

m 

O 

N 





m 

m 

m 

tn 





CM 

» 

« 

♦ 

# 

9 




E 

m 

CO 

CM 


m 



m 

S 

CM 

m 

f0 

m 

m 



« 


tM 

iH 

HI 

H 













1 











CD' 

ffl 

CM 

III' 






© 

O 

o 







# 

• * 

« 

• 





I 

CD 

m 

f^ 

IN 





i 

OD 


Ci 

ml 






m 

01 

CM 

o 

0\ 

H 







m 

m 

m 



CM 

9 

* 

« 

« 

i 

♦ 


o 

E 

t> 

ci. 

CM^ 

m 

10 



# 

s 

CM 


m 

; ' |f| 

10. 

10 


i 



rC 

Hi 

H 

H 

H 


% 


CM 









O 

10^ 







1 

i 

« 







8 

o 

CH 








CD 

r- 








Ml* 

o 


CM 

m 

10 





CM, 

If ■ 

HI 

m 

10 



-c4-*> 

i » 

* 

i 

♦ 

f 

# 


ufl 

H 

' o 

O 

CM 


Ml* 

10 


# 

E 

! m 

O 

' m 

OD 

in 

in 


O 

£ 


m 

10 

10 


#0 

C0 


% 


III' 

CM 


DD 

■ H 


o 


1 ^ 


r- 

m 

OD 

f0' 


cys 


! 

$ 

t 

# 

» 

' ♦ 


#. 


? m t 

Id 

10 

H 


O 


o 


1 

CD 

1# 

■fl 

m 

■ m 


<0 


1 

fMi' 

9f4 

H 

94 

H 


rl 


and 






CHAPTER ~ IV 


ANALYSIS OF COOPLED LAMINATES 


In this chapter a study has been made to analyse 
coupled antisymmetric angle ply laminates. Using success- 
ive iteration technique, A detailed investigation has been 
carried out with vario-us computational models for xmiaxial 
buckling problems with x edges simply supported and y 
edges clamped. 


4.1 Q0¥ERI^II®3 EQUATIONS 

The eqn. (2.15) can be rearranged for uniaxial 
bwkling problem as. 
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It may be ramentioned that 

e^# e2 »**.»*,e^g are defined in eqn* (2#8d) 


N 


with 


K =* w 0 


Prime denotes the differentiation with respect to *| • 

t 

Integrating eqn. (4.1a) successively one gets# 
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Integration of 4qn. (4 .lb) and eqp. C^i.lc) yielda 
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The constants of integrations k 2 _* kg# k^* kg# kg# 
k 7 # kg are to be determined from the boundary conditions. 

4.2 RAYbEIGH’S QUOTIIKT 

Using e<ja. (2.18) and e^* (2.20a)# the Rayleigh's 
quotient oan be written as 

X =• 

where 
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4.3 BOUNDARY CONDITION 

Using the relations {2.14a) to (2.140) in eqn. {2.13c) » 
boundary conditions for clamped y edges are obtained as 


u{0) 

u(l) 


w y (0) 
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m ■ W ( 0) 

m w{l) 


» w * (0) 0 

w w* (1) » 0 


{4.4) 


The constants of integration can be obtained using (4.4) in 
(4.2b) to (4. 2h). Defining the definite integrals as# 
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Constants of integrations are 


*« 12r2 - 6rj_r k2 = 21:^ - 6r^J k^ » Oi k^ * 0; kg • 


kg * Oi Itj m -r^i kg » 0 


C4«Se) 


4 . 4 . GEHBKAL STEI® IN SyCCESSX?E lUiaAiriON TKCHHIQ'Q® 

1, For the given fibre and matrix properties E^# Eg# 

and Y 2^2 generate [«J and fcj matrices nsir^ eqn»C2*4) 
and ^..Ceqn# (2*5) respectively* 



For the chosen ■© and thicknesses of the laminate# 
compute the coefficients of the stiffness matrices 
^ij# and D^j from eqns»C^*^^' 

For the specified aspect ratio (p)# length to thick- 
ness ratio Cr) and given mode ntmber# calc Til ate 
the constants g^# Z 2 # C^# c^# Cg# Cg, c^# Cg# 

®10* ®‘5P^* (4. id)* 

Select some Initial ftanctions w,v#u (say wO^) = v(»]) 
m u(y]) « 1) Which may or may not satisfy the 

boundary conditions and an arbitraary eigen value 
(say X » 10)* 

Find the value of eg from e<pi, (4. Id)* 

The constants of integration kg# k^# kg# kg are zero* 
Evaluate the constants k^# kg# kg and kj from the 
eqns* (4,6a) to (4*6e)* 

Using the sane already obtained or initialised 
function, wj. K* V u^ and Ujj# independently 

evaluate wj+it Vl* ^n+1' ’^n+l' ^tl '^ntl 

fr«M» eqns, (4* 2b) to C4,2h)* The suffix n denotes 

the iteration number* 

normalize the functions w«^^# v^^^# v^^^^# 

'^+1* \tl usiE^ a scalar factor (say 

Determine the eigen value X from the eqp, (4.3a)* 

Return the control to step S and terminate the iteration# 

vhen I ""^n^'^^n ^ ~ 6' 

(say ^ * 0*001)* 


a small number 
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11# For various values of mode ntrnbers# find tlie eigen 

val UBS The minimum eigen value out of these 

eigen valued is the buckling load* 

4.4*1 NUMERICAL STIDY 

Numerical computation was carried out on QFI^ six 
layered regular (equal thicknesses) laminate, with same 
material properties as mentioned in Chapter III, Table 4*1 
shows the eigen values obtained in the iteration process 
for m * 1 and 2 and buckling loads for various fibre orien- 
tations ®, A mesh size of 0*01 was selected for quadrature 
using trapezoidal rule. Following observations can be made 
from table 4*1, 

1* For ® * 20 and m « 1, the Iteration converged at an 

eigen value 97*91, For m = 2, it converged at 171*60* 
The buckling load 97.91 has an error of about 78% 
compared to Sharma*s ^8^ • 

2. For O » 30 and m » 1, the iteration converged at 

an eigen value 110.78* For m * 2 it converged at 
1S2.72. The buckling load 110.78 has an error of 
about 73% compared to Sharraa*s [l8“] 64 . 

3. For 9 m 46® and ra » the iteration converged at 
an eigen value 131.72, For m * 2 it converged at 
142.63* The budding load 133,72 has an error of , 
about 19*53% compared to Sharma's |_1^ 96* 
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Load. iiKiireases as © incirsases for 45^» This 
general trend agrees with the correct resiilt 

5. As the error involved is very large a detailed 

investigation in the case of © « 20° was carried 
out. 

Though convergence in eigen value was obtained in 
all the above cases tested, the value of the buckling loads 
are not correct. Various tests had to be made to check 
this behaviour.jj A test could be made to see the convergence 
behaviour of the iteration process# by using different starting 
functions say w(r)) m sin^r)} u(q) =» SinTty)#* v(v^) « SinTrrj 
which satisfy the boundary conditions? ii) A test cotild be 
made to see the possibility of roxand off errors that may 
occur in the iteration process*, by using double precision 
facility throughout the program? iii) A test could be' made, 
to see the influence of starting eigen value, by xising 
the eigen value X obtained from Rayleigh's quotient eqns. 

C4«|a) to (4 •31) for the initial functions w(r^) = v{r^) m nirj =» 

as starting value rather than previously assigned 10) valuei 

iv) A test eoxild be made to see the influence of the mesh size, 
by selecting laesh interval# say 0*001 instead of 0*01 in 
the quadrature using trapezoidal rule*: ' 

Table 4.2 gives the eigen values for the following 

casest 

i) convergence of iteration was achieved even xdien 

different starting fuixjtion were used. The convergi^ 



eigen value 0 \^) for the starting functions 
w(r|) » v(y|)= 1 was 97.91 and for 

w =a sin^vj? V a SinKif|i u =* SinT^i'^ the 
converged eigen value (X^) was 97.79, thoxagh 
the initially obtained eigen values are differ- 
ent from the previous case. This indicates that 
the iteration process is moving smoothly and is 
independent of starting functions. 

Computations in doxible precision do not affect 
the res\iits« ©le large discrepancy in results 
does not seem to be occurring due to round off 
errors. (X 3 ) values converges to 97*91. 

For a starting eigen value ^ « 50, the obtained, 
converged eigen value was 97.91 and for the 
second case using the starting eigen value got farom 
the Raleigh's quotiaat equation was {X^) 97.58* 

This Indicates that the Iteration process is 
Indei^ndent of starting eigen val ue and it is trying 
to achieve the same -convergence even tl^ugh the 
initially obtained eigen values in both the cases 
are different. 

The use of mesh size 0.001 instead of 0.01 gives 
the same iterated eigwi values and the converged 
eigen value (Xg) 97*91* 
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It may be presumed that evaluating the functions 

K+l* \+V K-i-V ^n+V ^^+1 '^n-f-l independently 

from eqns# (4 •2b) to {4*2h) usir^ the already obtained 

functions w^, v^# u^ and Ujj might have led to 

this incorrectness of the convergence# By using# the 
improved functions (say ^n+1^ addition to the 

already obtained functions 'Csay Wj^)# in the evaluation of 
functions (say from eqns. (4,2b) to C4,2h) one may 

expect better result. Various combinations are possible in 
the simultaneous iteration of fumtioas and various models 
could be studied, 

4,5 ITBRAIIOM SEQUENCES IN MCDELS 

The previous iteration scheme is denoted as model 1, 
In all the other models to be studied# only the steps 5 to 9 
of Model 1 had to be repl aced, 

4, 5,1 MOJEl. 2 

i) Choose the initial functions wCr)) = v(fj) = m(y)) » 1 
and an initial eigen value ^lO* 

ii) Evaluate u^^j^ end from eqns. <4.2g) and (4*2h)# 

using the functions wj# v^ and 

iii) Evaluate and from eqns* (4,2«) and (4,2f)# 

using the functions w^# w« a*^ \+l* 

Iv) Eealuate from eqns* (4.2b)# (4,2®) 

and C4*2d)# using the funetioas ^^ 4^1 

and Wjj. 



V) Normalise the function 

^n+l \+i the maxiram value of w , 

H’T 1 

vi) Find the eigen valued from eqn, (4*3 a). 


MC3DEI. 3 

Choose the initial functions u(y]) * vCr)) = u(l^) =1 
and an initial eigen value a= 10. 

Evaluate u^_j_j,and using the functions w«, 

v^ and u^. 

Evaluate v^^^^ and v^^^^ using the functions 
Normalise the furKstiona wJJ# w*# w , v * , » v 

*1 *1 H Hrrl. ^ l1iT*X 

'^n+l \i+l the maximum value of 

Find A from eqn. (4.3a). 

Evaluate '*'n+l the functions 

’^n+l* ^n+1 ^n* 

4.5,3 MODEL 4 

i) Choose the initial f unctions w (►) ) = v ('^) = u (rj ) «■ 1 

and an initial eigen value i, « 10, 

il) Evaluate w*^l* ^n+1 ^ntl functions 

V '^n 

lii) Evaluate and using the functions 

Vjj and Ujj. ■ 


4.5.2 

i) 

ii) 

ill) 

iv) 

v) 

Vi) 
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Evaluate and from eqns. (4.2g) and (4..2h), 

uains, the functions 

Normalise the functions w** a. w* m «• «■ 

n+i' ^n^l* ^n+1' V 

^n+1 ^ using the maxinmm value of w - • 

IUhl 

Find the eigen value \ » 


4 . 5.4 MODEL S 

^^®®se the initial fuiKstions w(r^) avCrj) au(v]) »1 
and an initial eigen value X = lo« 

ii) Evaluate wj+i* ^n+ 1 ' ^'n+l ^iug the functions u^, Ujj, 

V and w • 
a n 

ill) Evaluate u^^^ and using the functions 

“n+l' '"n “n- 

It) Evaluate v;;^^ and v^^j using the functions 

"n+l' % ’^,+1- 

V) Hoin«.liae ^ 

using maxirntmi value of 

vij Find the eigen value X from Rayleigh's quotient eqn. 

C4 . 3a} • 


4.5.5 MCDEL 6 

i} Choose the initial functions w{r]) » v^) « u(»^) » l, 

ii) Evaluate v^^j and v^^^ fro® eqps. (4.2e) and (4.2f), 

using the functions u^*„ 

ili) Normalise the functions w«, v^^l# v^^^* 

using the maximum value of w^j. 



Find the eigen valued • 

IJvaluate tising «ie fimctions 

^n* ^n-¥l '^n* 

Evaluate usir^j the functions 

^ntl' Vl «n- 

mJMERlSAI. STWY IN MODELS 

The same QFBP six layered regular 20° angle ply 
laminate was taken for numerical study in variotis models. 
Table 4*3 shows the eigen Talues obtained in the iteration 
process for the models* The starting function was selected 
as w(»^) ■« v(ri) * u(v^) « 1 . A mesh size of 0.01 was used 

in the quadrature using trapezoidal rule. 

following observations can be made from a study of 
the tablet 

i) In model 2, the Iterated eigen values oscillate 
between 64 to 78 in the initial five iteration and 
finally it converged simaothly at a value of 81.39. 

ii) In model 3# the iterated eigen values oscillated 
between 64 to 77 in the first five iterations and 
finally it converged at a value of 80.88.er 

iii) In model 4# the iterated eigen values oscillated 
between 69 to 86 in the first five iterations 

and finally converged at 79*58 after another three' 
iterations. 


iv) 
V } 

vi) 

4.6 
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ivj Im TOc!«l 5# the iterated eigan values oscillated 

between 70 to 84 in the first five iteratives and 
converged at 81*07 after another nine iterations* 

roodel 6# the iterated eigen valnes converged 
in five iterations to a value of 79.24* 

vi? Compared to model 1* all the above models are 

converging closely around 80* These models are 
better than model I, which converged at a value 
around 98. As these restilts differ by about 50% 
with the correct results, model 5 was further investi- 
gated* 


4.7 DETAILED HOMERlCAL ANALYSIS IN MODEL 5 

A detailed study could be made to analyze the eigen 
values in model 5 iterations scheme, by varying the nun^er 
of layer# in OIRP 20® angle ply regular laminate* Table 4.4 
shows the eigen values obtained in the model 5 iterations 
process for varied ntradjer of layers from 2 to 50. 

i) In 2 layered laminate, the iteration converged at 
74.25. Sbaxma £l8^ obtained the buckling load as 
26, The error is about 185%, 

ii) For 4 layered laminate, the iteration converged at 
74,74* Sharma 

48, The error is about 55%, 


,8^ obtained the buckling load as 
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lit) 


v} 


vl) 


For 6 layerBd laninate# ■the iteration converged 
at 81. 07# Shaxma 18 obtained the bi^kling 
load as 55. The error is about 47%. 

For 20 layered laminate the iteration converged 
at 85*27. Sharma 18 obtained the buckling 
load for Infinity layers as 58. The error is 
about 46%. 

For 50 layered laminate the iteration converged at 

86 . 01 . 


As the number of layers are increased, the converged 
eigen value also increases gradually. This behaviour 
is logically true in laminate analysis. The coupling 
coefficients lj_g aiKi B 26 aero as the nwaber 

of layers are increased and the laminate, b^ave ■ 
as decoupled laminate. The error in 'percentage -is. 
reduced as the n’umber of layers' are increased in 
this model scheme, ; This indicates that the presence 
of and l 2 g terras has influence in obtaining the 
bucklii^ load. 

In thin classical plate theory, "the effect of thick*- 
ness to length ratio (r) is not felt. But the governing 
equations used in the present investigation appear to be 
dependent upon r. A study of the dependeaice of on r is 
pres«ated in Table 4.5. 
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Table 4,5 shows the (i) eigen values obtained for 
r « 1 case in 2o® six ply regular QFRP laminate with start- 
ing functions w(r)) a virj} m u(»^) a i, (u) r a 1 case with 
starting functions wCr^) a v(r^) a w(rj) a 1 and (iii) r » 2C 
case with starting functions w(k|) a Sin^i^ I v^q) « sinfCi^ ; 
u(<^) m sin7TJ|, The following observations could be made: 

value of r is not at all influencing the conver- 
gence of the iteration process. For r * 1 , and 20 
and for different starting functions the eigen value 
converged to same 81. 07, 

Table 4.6 shows the values of u, v and'w displace- 
ments at various discrete points and the constants of inte- 
gration 1 ^ 2 , ki and for the final two stages of iter- 
ation process in six layered regular tSFSP 20®' laminate, 

Thoij^h the w valties converged at the final two stages 
of iteration, the u and v values did not converge. The 
constants of integration k^* ^7 also not 

converging because of the non-convergence of u, v. The 
incorrect buckling load might have obtained due to this 
behaviour. 

Table 4*7 shows the values of the particular integral 
terras H 2 # 3^* ®^ 4*****%1 appearing in the kayl eights 
quotient e< 5 pa* (4* 3 a) for the final two stages of iteration 
in six layered"(|?fUP 20® laminate. - Convergence in each terras 
Hj,# H 2 # Hj# 1 %* H 7 # Hq, and are achieved 

for final successive stages <Sfi iteration as observed in 
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table 4«7» This indicates tha^- 

^nat GOiwergence occurs in each 

individual particular integral 4 . ^ . 

^rai terms of the Rayleigh’s 

quoti ait expression in additi nr. * , 

^ion to the eigen value • 

The solution teehnim.« 

*3^6# appears to give the solution 

of the eqns. {4.1a) to (4-.lei „ 

• Howeveri, the iteration conver- 
ges to incorrect results. no4.^i, ^ ^ 

^©tailed investigation reveals the 

convergence on the eigen velue, bet the mode ehape, epeclelly 
Inplane dtepleoemente n and y, eomrerge properly. It la 

difficult to pinpoint the reasona. Probably this is due to 
inherent limitation of the method. One of the reasons could 

be seeking of the convergenoe on three variables u. v and w 

through iteration* 


It is possible to reduce the present set of governing 

equations into s set of two ordinary differential 

equation in v and w, by eliminating the variable u. The 
solution then boUs down to .elution involving only two 
unknown dependent variables, mis approach la presented in 
the next section. 


4.8 SOLUTIOHT USIJJG TWO mRlAl|,gg r and w 

The governing equations Ci.la) to (4.1c) are rewrittw 
here for clarity sake. 

w'- » CjV + c^u» + CgV- + Cg U-' (4.1a) 

V- « ♦ «8^ + CgV ♦ C3L0«’ <4.1b) 

Proa e<|ii* (4.1b)» tme gets* 


(4.1c) 


Klimlnation of u" from eqns, (4,1c) and (4,6) gives 

n m v»»« + ^ fjW"* + f^w» (4.7) 

+ f^v- + + f^w- (4.8) 

Substituting u‘ from egn. (4.4) into eqn. (4. 2}» one gets 
IV TV 

V * gj^v" + gjV + + g^ w« ' + ggW** + g^w 


(4.9) 

It may be noted that in the above equation, following constants 
have been defined! 
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f. 


^6 


^4 

IJ' 


Integrating egn. (4.9) successively# one obtains 
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V" » V + I I vd if + w” *1* g^w* + ggW + gg ^ f wd^ 

o o 

+ + k^Q ^ 4 ^i 2 ) 

-■ - ai I rd^ . g, + 93 «■ * 94 W 4 95 P w dT 

w \d# v o 


V « 


i o ®2 /f /o''‘^'l* ^3“ * 34 P'^ 

+ 95 i* /w a>)^ + gg ,1 x'i/l_/Iw a^‘‘ + kg §- 

00 o o o o 


+ k 


10 


■¥ 


+ k 


11*1 


(4.14) 


listing th« boundaty conditions (v(l) »» n(o) « nCl) *» 0)# 
constants kg, k^g# can be evaluated as a solution of a 
set of three linear algebric equations* Omitting the 
details the constants can be evaluated from the relations* 


kfs ** t 

Hq 


9 S*« ' ^0 " ®i* ^9 + 37? • *2 ■** ®3^9 


(4*15) 

where the following constants are required to determine kg# 
ki « and kt * I 

m» %0' 


m 




2 

1 

t 


m + 

g^ W** ( 0 ) 

(4.16) 

' ( 0)1 

% 

'jP' ' ''^' "' mi "'#1^ '' 4^s#" 

1 *2 

(4*17) 


i < . t g| 

w*** ( 1 ) ♦ ^ 4 . w" ( 1 ) 

1 

+ gg wd 

Q 

(4*18) 
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f vdtj + 
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o o i ■*■ ^5 

.nn_3 


36 { ll 


(4.19) 


*g “ 


+ ^2 ^®3 0*5} J 


M>f^ • in* 


(f^ »4 "*• ^2 ^^5 ^2^ ~*^f3 V"* (1)) 

f « 


(4,20) 


Sq » 


r1 trr^K^ 


^ n n 

t Al . 


i i '''^’' "■ “2 Xo i i i ''■^1 + % X ''^1 


+ 9$ -Iq I ■*■ 94 


r ^ 4 

i t (’ f »d^ (4.21) 


«9 “ -<*8 + »2 ■*■ 2^*' *10 " ®3 + + 


•8 


(4*22} 


It i» q^te evident that suteessive Integration on 
eqn* (4.9) can be carried ont t© d®tej^£jje 

We now need to eXinninate n froia the e^* (4*ia) for w, 
Using eqns* (4*lb) and (4«1©)» w rewritt«nt 

as 

+ ' Og^ * 4^-OggM* (4*23) 


where 


m d^ ir" + djV t djw" + 

^ m fg ir"* t f|ir» + fiW** 4, £ V* 


( 4 . 7 ) 

( 4 . 8 ) 


and constan**^, ■ 
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«1S ^-^11 ^4^5 e.) 

(sj + egCg) * 


®16 * ®7 - aSgA ? z„ « 


^^6^4 <=6^6glf^aiioC ) 

^®1 ®6®9^ 


* 


f 2 2. 

v at i-^ ) 


(^1 4- eg Cj) 


f c^ » 


(eg-t'cs.^e^ + ®7®g©|3^) 

(e^ + egCg) 


'18 


imj + egc^^ 4- c^Cq -i- egCgCj^^) 

^®1 + % 69 ) 


(4,24) 


«ie present problem, thns, boils down to the solution 
of fourth order eqna, (4,9) and (4*23) in v and, w, Por 
known v and w, n* and n which are required to, evaluate 
are determined from eqns. (4,7) and (4,8) respectively* 

4*9 NUMSRJCAh STTO W11?K tWO milftBIiES 

3oluticm was obtained for the laminate used in 
section 4,4*1 for O » 20 , Two models for iteration were ■ ■ 
used, . ■ ■ 

4.9.1 MODEL 1 

Iteration follows similar procedure as used in the 
case of three variables, h very brief outline is given 
below* 

1. Choose the initial funotions and X as 
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a* Detexroine r and then u from egn. (4,7). 

3, Using above computed v, u determine w and calculate X 

4, determine u and u*, 

5, Normalise u, v and w with respect to Go back 

to step 2 and continue iteration till converges. 

Numerical results obtained for m * 1 and 2 are 
presented in Table 4*8 which gives associated with each 
iteration cycle, ibr m 1, iteration does not converge. 
Eigen value finally appears to oscillate between 59,2 and 190 
(^PP^ox,) , However for m *» 2 convergence occurs in 3 iter-* 
ations* It was further observed that the mode shapes is 
displacement functions u» v and w converged accurately* 

4. 9,2 MODEL 2 

In the previous case v obtained after iteration was 
UB«wl in computing w. In this case w equation is integrated 
without using previously cal.culated v i,e» integration of the 
two equations is carried- out independently* For m » 1 
associated with first few iterations is presaited in Table 
4.8, Oscillations observed in model 1 were also seen in this 
case, 

4,10 aiNSRAL OBSENfATIOlli 

h relatively detailed investigatioa was made* The 
results Obtained are net presented here. However* following 
general comments based on the stud^ can be made* 



The value of X » 56.2 i.e. Model 2 ia the exact 
result. Value X » 59,2 is fairly dclose to this 

value. The iteration thus appears to oscillate 
hetM®m a correct and highly incorrect value. 

Correct value was obtained within first few iter- 
ations. One tends to believe that the method yields 
reasonably correct results in 4 to 5 iterations, SO 
using model 2, a study was made for' different values 
of ©• Results do not substantiate the conclusions 
reached for © » 20®. 

For m « 2# precise convergence occurs for © m 20®. 

It is a highly pmplexing behaviour as no other ^ 
change# except in -the value of m# has been made. 
Probably the solution process is very sensitive to 
the values of the coefficients of the derivatives 
on the right hand side of the equations (4.9) and 
(4*23). Oscillations were c^served for m « 2 when - 
© was varied. ' ' 

It appears that siiscessive iteration in such problems 
may lead to correct result. It has# however, failed 
to give the desired solution in the present study. 
Even if the method may worh in some cases# its 
reliability when the nmBmr of dependent variables is 
more than 1 appears to be doubtful, fh# method# 
therefore# should be att«ipted when the dependent 



variable only U 

I*®Sicaiiy, method seems to be acceptdsle* Inves- 
tigation attempted In the present study are unable to 

pinpoint the reasons behind the failures of the 
method. 
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concLnsioiia 


The preaent work Is e moa«:ate contribution to the 
analysis of ooo(posit« plates. Based on the results obtained 

certain coaclteicua cf tn be drawn. 


la* 


b. 


c. 



Succeasiire iteration technique works well in the 
buckling and free iribration probleitis of decoupled 
thin antisyw^etric angle ply laminates. The 
technique giwes slightly higher elgwi yalues ' 


contpared to other solution techniques* 

The coa^putetiomel efficiency is high and the solution 
cowerges closely with the results of other technique** 
la meiy of tlwi eases ^*e soli^icm is obtained in atmost 
fiye iterati«is* 

Prograsming requires less effort and for differ«mt 
boundary conditions the changes to be made are only 
in finding the waliws the oonstants of integration. 


A mesh cite of 0.01 InteurdL in quadrature using 
trapezoidal rule is sufficient* Comrergeace. cirlteria 


^ ^n+1 "* \ ^^n-flll ^ the two omcessiye 

iterations to-toiwerge .in sutficieat for aocuracy. 

^h® techniqttie «oaire«fesi;iKi>|ii^,s«ie:.eig,«yalue* tren.if 
differmit inAtlal 



In thin conqposite plate hi* 2 kling problems with sa, sC 
«nd CC boundary conditions^ optimization can be made 
with respect to fibre orientaUon using this technique. 

To imjrease the load carrying capacity in QWt aw and 
^ conditior^ the laminate has to be designed for # m Q 
fibre orientation. The buckling load increases as the 
aspect ratio ii«reases* For various biaxial loading 
ratios interaction curves can be drawn with this technique* 
which will be useful in design* 

In free vibration problems with CC condition the plate 
has to be designed for 0 a go® and in FW condition for 
• » 0 to increase the natural frequency* hs the aspect 
ratio ijxsreases# the natural frequency decreases* reaches 
minimum and then increases in GC and SS condition whereas 
it drops continuously in FF condition. 

The technique can be applied readily* In optimal design 
of laminates* as it consumes less computational time. 

In optimisation process the analysis has to be repeated 
many times for various parameters* 

Optimization studies can be made in hybrid thin anti- 
symmetric angle ply laminates* as it is possible to 
analyse decoupli^ thin laminates* using this technique. 

Ttm present study has been made* using small deformation 
theory* Optimisation stiidies may be performed in lazge 
deformation domain* 
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Fibra TOliime fractions in addition to fibre orientations 
of laminae catn be taken in optimal design. 

2* Successiire iteration techniquet when attempted in thiio. 

coupled antisymmetric laminate buckling problem failed 
to give correct eigen value X , oven though positive 
convergeiwe in X io achieved. Various model schemes have 
been tried to investigate the failure of the method. 

The method se«ns to be doubtful when the number of depen- 
dent variables more than one appear in the governing 
equations. The solution process may be sensitive to 
the values of the coefficients of the derivatives^ appear- 
ing in the governing equations* The present study is 
unable to pinpoint the reasons behind the failure of 
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